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Abstract 



The one-dimensional Navier-Stokes equations are used to derive analytical ex- 
pressions for the relation between pressure and volumetric flow rate in capillaries of 
five different converging-diverging axisymmetric geometries for Newtonian fluids. 
The results are compared to previously derived expressions for the same geometries 
using the lubrication approximation. The results of the one-dimensional Navier- 
Stokes are identical to those obtained from the lubrication approximation within a 
non-dimensional numerical factor. 
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1 Introduction 

Modeling the flow through capillaries of converging-diverging geometries is an im- 
portant subject and has many scientific and industrial applications. Moreover, it 
is required for modeling viscoelasticity, yield-stress and the flow of Newtonian and 
non-Newtonian fluids through porous media [1-6]. 

There are many previous attempts to model the flow through capillaries of 
various geometries. However, they either apply to tubes of regular cross sections [7, 
8] or deal with very special cases. Most these studies use numerical mesh techniques 
such as finite difference and spectral methods to obtain numerical results. Some 
examples of these attempts are Kozicki et al. [9], Miller [10], Oka [11], Williams and 
Javadpour [12], Phan-Thien et al. [13, 14], Lahbabi and Chang [15], Burdette et al. 
[16], Pilitsis et al. [17, 18], James et al. [19], Talwar and Khomami [20], Koshiba et 
al. [21], Masuleh and Phillips [22], and Davidson et al. [23]. 

In this article we use the one-dimensional Navier-Stokes equations to derive 
analytical expressions for the flow of Newtonian fluids in tubes of five axisym- 
metric converging-diverging geometries, some of which are schematically depicted 
in Figure 1, and compare our results to previously derived expressions using the 
lubrication approximation [24]. 



Figure 1: Profiles of converging-diverging axisymmetric capillaries. 
The widely-used one-dimensional form of the Navier-Stokes equations to de- 




1 INTRODUCTION 



6 



scribe the flow in a tube of length L where its axis is aligned with the x axis and its 
midpoint is at x = is given by the following continuity and momentum balance 
relations respectively assuming negligible gravitational body forces [25-37] 



dQ d 
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In these equations, A is the tube cross sectional area, t is the time, Q (= Au) 



is the volumetric flow rate, x is the axial coordinate along the tube, a( 



Ju 2 dA, 
Au' 2 > 



IS 



the momentum flux correction factor, p is the fluid mass density, p is the pressure, 
and k is a viscosity friction coefficient which is given by k = with p being 

the fluid dynamic viscosity. For steady flow, the time terms are zero, and hence Q 
as a function of x is constant. The momentum equation then becomes 
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For a flow in the positive x direction, the pressure gradient is negative and 
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This expression is dimensionally consistent. 




Figure 2: Schematic representation of the radius of a conically shaped converging- 
diverging capillary as a function of the distance along the tube axis. 
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For a tube of a conical profile, depicted in Figure 2, the radius r as a function 
of the axial distance x is given by 



r(x) = a + b\x\ - L/2<x <L/2 



(12) 



where 
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Hence, Equation 11 becomes 
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1.2 Parabolic Tube 



For a tube of parabolic profile, depicted in Figure 3, the radius is given by 



r(x) = a + bx 2 - L/2 < x < L/2 



(19) 



1.2 Parabolic Tube 
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Figure 3: Schematic representation of the radius of a converging- diverging capillary 
with a parabolic profile as a function of the distance along the tube axis. 
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a = R r . 
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Therefore, Equation 11 becomes 
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1.3 Hyperbolic Tube 

For a tube of hyperbolic profile, similar to the profile in Figure 3, the radius is 
given by 



r(x) = Va + bx 2 - L/2<x<L/2 a,b>0 



(24) 



where 
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1.4 Hyperbolic Cosine Tube 

For a tube of hyperbolic cosine profile, similar to the profile in Figure 3, the radius 
is given by 



r(x) = acosh(6x) — L/2 < x < L/2 



(29) 



where 
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Rmin and b = — arccosh ( ™ ax ) (30) 
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1.5 Sinusoidal Tube 




Figure 4: Schematic representation of the radius of a converging-diverging capillary 
with a sinusoidal profile as a function of the distance along the tube axis. 



For a tube of sinusoidal profile, depicted in Figure 4, where the tube length L 
spans one complete wavelength, the radius is given by 



r(x) = a — b cos (kx) 



- L/2 < x < L/2 



a > b > 



(34) 
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where 
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Hence, Equation 11 becomes 
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On performing this integration, the following relation is obtained 
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On taking lim L - I the following expression is obtained 
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Since 5 < —1, p > as it should be. On substituting for B, b and in the last 
expression we obtain 
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It is noteworthy that all these relations (i.e. Equations 18, 23, 28, 33 and 
42), are dimensionally consistent. Moreover, they are identical to the lubrication 
approximation formula (refer to Table 1) for a = 4/3. 



Table 1: Lubrication approximation table. These formulae are derived in [24]. 
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2 Conclusions 

In this paper we derived analytical expressions relating the pressure drop to the vol- 
umetric flow rate for Newtonian fluids in five different converging-diverging geome- 
tries using the one-dimensional Navier-Stokes flow equations in axisymmetric tubes. 
The results obtained in this paper are identical, within a non-dimensional numer- 
ical factor, to those derived in [24] using the lubrication approximation. These 
expressions can be used in various practical scientific and engineering situations 
to describe isothermal, uniform, laminar flow of incompressible, time-independent 
Newtonian fluids. These situations include the flow in corrugated vessels and the 
flow in the pores and throats of porous media where the converging-diverging na- 
ture can be idealized by one of these simple geometries. The analytical method 
can also be used to derive expressions for axisymmetric geometries other than those 
presented in this paper. 
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Nomenclature 

a correction factor for axial momentum flux 

k viscosity friction coefficient (m 2 .s _1 ) 

H fluid dynamic viscosity (Pa.s) 

p fluid mass density (kg.m~ 3 ) 

A tube cross sectional area (m 2 ) 

L tube length (m) 

p pressure (Pa) 

Q volumetric flow rate (m 3 .s _1 ) 

r tube radius (m) 

Rmax maximum radius of converging-diverging tube (m) 

Rmin minimum radius of converging- diverging tube (m) 

t time (s) 

u local axial fluid speed (m.s -1 ) 

u mean axial fluid speed (m.s -1 ) 

x axial coordinate (m) 
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